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■ Abstract 

• We analyse the linear kinetic transport equation with a BGK relaxation operator. We study the large 

^ ■ scale hyperbolic limit {t,x) {t/e,x/e). We derive a new type of limiting Hamilton-Jacobi equation, 

pL^ , which is analogous to the classical eikonal equation derived from the heat equation with small diffusivity. 

We prove well-posedness of the phase problem and convergence towards the viscosity solution of the 
Hamilton-Jacobi equation. This is a preliminary work before analysing the propagation of reaction 
fronts in kinetic equations. 

' Resume 

. ■ Une equation eikonale cinetique. Nous analysons une equation cinetique lineaire de transport avec un 

' operateur de relaxation BGK. Nous etudions la limite hyperbolique de grande echelle {t, x) — >■ {t/e,x/e). 

, Nous obtenons a la limite une nouvelle equation de Ifamilton-Jacobi, qui est I'analogue de I'equation 

eikonale classique obtenue a partir de I'equation de la chaleur avec petite diffusion. Nous demontrons 
le caractere bien pose de I'equation verifiee par la phase, ainsi que la convergence vers une solution 
de viscosite de I'equation de Ifamilton-Jacobi. Ceci est un travail preliminaire en vue d'analyser la 
propagation de fronts de reaction pour des equations cinetiques. 

>: 

(N 

^! Version frangaise abregee 
(N. 

, Nous considerons un modele cinetique lineaire avec un operateur de relaxation BGK, pose sur un ensemble 

(<~^ ■ de vitesses V symetrique et borne. On analyse le comportement de I'equation dans la limite hyperbolique de 

^ I grande echelle (t, x) (|, |), 



>< 



dtf + v-V^r ^-{M{v)p' ^ n, {t,x,v)eR+xWrxV. (1) 
e 



$H I Nous demontrons que la phase ip^ definie par la relation f^(t,x,v) = M{v)e * ; ' converge (localement) 



uniformement, lorsque £ — >■ 0, vers une fonction ip'^{t,x) independante de v. De surcroit, la fonction ip^ est 
solution de viscosite de I'equation de Hamilton-Jacobi suivante, 

^^''^ dv = l, (t,x)eR+xM". (2) 



V l-dtip°it,x)-vV^ip°it,x) 



Cette equation pent se reecrire sous la forme canonique dtp'^ + H{y x^^) — pour un hamiltonien effcctif 
Hij)) qui est lipschitzien et convexe. En regie generale, nos travaux consistent a homogeneiser I'equation ^ 
par rapport a la variable de vitesse. Le probleme aux valeurs propres dans la cellule V s'ecrit comme suit : 
trouver un vecteur propre Q{v) et une valeur propre H(j)) tels que 

(1 + H{p) -v-p) Q{v) = / M{v')Q{v')dv' . 
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La demonstration du passage a la limite de ((TJ vers ^ s'appuie sur une serie d'estimations a priori qui 
demontre que (f^ appartient a I'espace de Sobolev W^'°°, avec un controle uniforme en e > (Proposition 
l2.1l ci-dessousV Dans un deuxieme temps, nous demontrons que toute fonction test ifj'^{t,x) de classe telle 
que (p^ — ijP admet un maximum local en {iP ,x^^, verifie 

M(v) 

^ ' -dv <1. 



Ceci demontre que (p^ est une sous-solution de viscosite I'equation de Hamilton- Jacobi ([2]). Un raisonnement 
identique montre qu'il s'agit aussi d'une sur-solution de viscosite. La demonstration se base sur la construction 
d'un correcteur microscopique r]{t, x, v) defini de fagon ad-hoc par la relation 



1 Large-scale limit and derivation of the Hamilton-Jacobi equa- 
tion 

We consider the following kinetic equation with BGK relaxation operator: 

dtf + v-V^f^M{v)p~f, it,x,v)eR+xW'xV, (3) 

where f(t,x,v) denotes the density of particles moving with speed v Cz V at time t and position x. The 
function p{t, x) denotes the macroscopic density of particules: 

pit,x)^ / f{t,x,v)dv, (i, a;) e R+ X M" . 
Jv 

Here V denotes a bounded symmetric subset of M". We assume that the Maxwellian M is symmetric and 
satisfies the following moment identities: 



M{v)dv^l, / vM{v)dv^O, / v^M{v)dv 



V 

In this paper we focus on the large scale hyperbolic limit (t, a;) — > (|, |), e — ?► 0. The kinetic equation (jSj 
reads as follows in the new scaling: 

dtf + v-^^r- ^-{M{v)p' - n, {t,x,v)en+xW xV . (4) 

£ 

Clearly, the velocity distribution relaxes rapidly towards the Maxwellian distribution. This motivates the 
introduction of the following Hopf-Cole transformation: 

r{t,x,v) = M{v)e-'^^^ . 

where we expect the phase to become independent of w as £ — > 0. To avoid technical complications due 
to ill-prepared data, we set Lp^{Q,x,v) — Lpo{x) > as an initial data for The equation satisfied by Lp"^ 
reads: 

dtip'' + v - V^p^ ^ j M(u') (l -e'^^^) dv', (t, a;, w) G M+ x M" x F , (5) 

Theorem 1.1 Let V C M" he hounded and symmetric, and M e L^{y) he nonnegative and symmetric. 
Then ip^ converges (locally) uniformly towards p^ , where p'^ does not depend on v. Moreover pP is the 
viscosity solution of the following Hamilton- J acohi equation: 

^^^"'^ dv = l, (t, x)£R+x R" . (6) 



Iv l-dtp°it,x)-vV^p°{t,x) 
The denominator of the integrand is positive for all v Cz V . 
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The last observation in Theorem ll.il is not compatible with an unbounded velocity set. 
One can understand this as an homogenization problem in the velocity variable. Moreover, the associated 
eigenproblem in the cell V writes: Find an eigenvector Q{v) such that: 

(1 - dtip° - V ■ V,^") Q{v) = / M{v')Q{v')dv' . 



V 



This eigenproblem can be solved explicitly, and yields formula ([5]). 

Thanks to monotonicity properties, we can boil down to the classical framework of first order Hamilton- 
Jacobi equations. Indeed, writing equation ([5]) as G{dt^p^ ,V x'P^) — 0, we observe that G is increasing with 
respect to the first variable. Hence the equation is equivalent to 9t(p° + H(y x^^) — 0, where the effective 
Hamiltonian H is defined through the implicit formula. 



Differentiating ([7]) we obtain. 



V il + H{p)-vp) 
M{v) 



iVH{p) -v) dv^O. 



We deduce ||V-ff||oo < Knax- This is in accordance with the underlying kinetic equation, since VH can be 
interpreted as the group speed, which is bounded by the maximal speed of the particles. Differentiating ([7|) 
twice we obtain 

^"^^""^ ■ dv^ B^H{p) = 2 / JfP- ^ (VH(p) -v)<S) {VH{p) ~ v) dv . 



y (1 + H{p) ~vpf J Jy (1 + H{p) - V ■ pf 

We deduce that the effective Hamiltonian is convex. 

As an example, we can compute the effective Hamiltonian H in one dimension for a constant Maxwellian 
M=\onV^ (-1, 1). We obtain H{p) = It is equivalent to 6^\p\'^ for small p {9^ = i). Another 

example where the effective hamitonian is explicit is given by the Maxwellian M{v) — ^ {Si + S-i), though 
it is not a function. This corresponds to a two velocities model (also known as the telegraph equation, see 

In this case we obtain the relativistic hamiltonian H{p) = V^+^p 1^ 
Interestingly enough, we obtain a Hamilton- Jacobi equation which differs from the classical eikonal equa- 
tion. The latter could have been expected from the following argumentation. The formal limit of equation 
(jl]) at order 0(e) is the heat equation with small diffusivity: 

dtp' = sd^Axp' , [t, x) e M+ X M" . 

It is well-known that the phase (p'^ = —e log p' satisfies in the limit e — )■ the classical eikonal equation in 
the sense of viscosity solutions [5l [TOl El [71 [11] : 

5t0° + 02|V,0"P=O. (8) 

We only have asymptotic equivalence between the two approaches for small \p\ as can be seen directly on 
by Taylor expansion: H{p) ~ 

In Figure[l]we show numerical simulations of the kinetic eikonal equation ([6]), with a constant Maxwellian 
on y = (—1, 1), and we compare it with the classical eikonal equation ([8|). 

We end this introduction by listing some possible extensions of Theorem 11.11 for other choices of transport 
and scattering operators. We will develop a more general framework in a future work. 

1. In the case of the Vlasov-Fokker-Planck equation, 

dtf + V ■ Vxf - VxV{x) ■ V^f = -V, ■ + vf) , (t, X, v)eR+x M" x M" , 

e 

we obtain simply the eikonal equation dt(f>'^ + IV^^^^P = in the WKB expansion — Af(u)e~^, 
where M{v) is a Gaussian. 
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Figure 1: Numerical simulations of the Hamilton-Jacobi equation 9f(p + H{dxf) = 0. {left) The kinetic 
eikonal equation ^ where M(v) = i). (right) The classical eikonal equation H{p) ~ 6'^\p\'^ ([S]). In 

both cases the initial condition is a parabola. 



2. It is challenging to replace the BGK operator in Q by a convolution operator L{f) — K * f — f, where 
iC is a probability kernel [T] . However in this case we are not able to solve explicitly the eigenproblem 
in the cell V. 

3. In a forthcoming work we will investigate the propagation of reaction fronts in kinetic equations, 
following [7l[9lll[2]. 

2 Proof of Theorem 11.11 

First let us mention that the solution ip"^ remains nonnegative for all times. We proceed in two steps. First 
we prove uniform estimates with respect to e > 0. It allows to extract a uniformly converging subsequence. 
Second we identify the limit as the viscosity solution of equation © using the maximum principle. The 
second step relies on the construction of a suitable corrector ri{t,x,v) [61[8]. 

Step 1. Existence and uniform bounds. 

Proposition 2.1 Let V C M" be a bounded subset. Assume M e L'^{V) and ipo G iyi'°°(]R"). The 
kinetic equation ^ has a unique solution (p^ € W^'°° (K+ x M" x V). Furthermore, the solution satisfies 
the following uniform estimates: 



o<^'{t,-)<yo\\oo, (9) 

||V:,¥''(i,-)l|oo < ||V^(^0|loo, (10) 

||V„(^'^(t,-)||oo <i||V:,(^0||oo, (11) 

\\dt(p''{t,-)\\oo < Kiaxll (12) 



Proof. We obtain a unique solution ip^ from a fixed point method on the Duhamel formulation of ([5]): 

^%t,x,v)^ipQ{x-tv)+ J M{v') (l ~ 6 "'"'°'°"°"'"'"/""°'"'°"'"'' j (13) 
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We obtain directly, 

Ve > 0, 0<ip''{t,x,v) <ipo{x-tv) + t. 

This ensures that ip'^ is uniformly bounded on [0, T] x M" x V. To prove the bound ©, we define ipg{t, x, v) = 
ip'^{t,x,v) — 6t — (5^|a;p. For any 6 > 0, ipg attains a maximum at point {ts,xs,vs). Suppose that tg > 0. 
Then, we have 

dt(p%ts,xs,vs) > 6, V xip^{ts-,xs,vs) ^25^xs. 
As a consequence, we have at the maximum point {ts,xs,vs): 

0>J M{v') ^1 - g ^-'"^-"^-"^';'^^"^-"-'-" >s + 2vs5'^xs > 5 ~ 2V^^j'^\xi, \ . (14) 

Moreover, the maximal property of {tj^xs^v^) also implies 

W\\oo~5^\x5\'' > v'{ts,xs,vs)-Sts-d^\xs\'' > f'{0,0,vs) > 0. 

We obtain a contradiction with ([T4|) since / {2Vinayi) < \xs\'^ < i5^*||<(2'^||oo cannot hold for sufficiently 
small (5 > 0. As a consequence ts — 0, and we have, 

V(i,a;,w) G [0,T] x M" x V, ^'{t,x,v) < ip° {xs,vs) + 6t + 6^\x\^ < \\ipQ\\oc, + 5t + 6^\x\^ . 

Passing to the limit (5 — > 0, we obtain (|9|). To find the bound ([T0|) . we use the same ideas on the difference 

(pf^{t,x,v) = ip^{t,x + h,v) — ip^{t,x,v). The equation for tpf^ reads as follows, 

dt^l + V ■ V,vp| = ^ M{v')e^^ (^1 - e^^^ dv' . 

Using the same argument as above with a (5— correction, we conclude that 

y{t,x,v) e[0,T]xR" xV, ipl{t,x,v)< sup \ip°{x + h,v) - ip°{x,v)\ 

{x,v)eMxV 

The same argument applies to ~(p%, 

dt i-^l) +vS/, i-^l) = - M{v')e^ (^1 - e-^^^^^^M^ dv'. 

so that the r.h.s has the right sign when —(p% attains a maximum. Finally, 

\f{t,x,v)e[0,T]xR'^ xV, \^i{t,x,v)\< sup l^^^ix + h,v) - ip"{x,v)\ < \h\. 

(x,v)£RxV 

from which the estimate (|10p follows. 

To obtain regularity in the velocity variable (ITT]) , we differentiate ([S|) with respect to v, 

{dt+v- Vx) iV^ip') = -g,{ip'')V^ip' - V^(p^ 
where gei^fi^) — ^ Jv M{v')e^ / dv' > 0. Multiplying by ; we obtain 

{dt+v V.) (|V„^^|) = -<?,(^^)|V„¥.^| ~ (^V,^^ • ^^yj (15) 

< llVaj-^olloo ■ 

from which we deduce (jlip since V^^^o = by hypothesis. 

Finally, the bound (|12p is obtained similarly as the bound on VxV?^ (fTO|). using the difference x, w) = 

if' {t + s^x^v) — ^p^ {t,x,v). We obtain 

V(t, X, w) e [0, T] x M" x y, |(^^(t,a;,w)| < sup |(/3^(s, x, w) - -^"(x, w)| . 

(x,t))eRxy 
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We use the Duhamcl formulation to estimate the last contribution: 

\(p''{s,x,v) - ip°{x,v)\ < \(po{x - sv) - (poix)\ + o{s) . 

The estimate follows. □ 
Step 2. Viscosity solution procedure. 

From Proposition l2.1l we deduce that the familly (f^)^ is locally uniformly bounded in W^'^ (R+ x M" x V). 
Then, from the Ascoli-Arzela theorem, we can extract a locally uniformly converging subsequence. We denote 

by the limit. Furthermore, from the fact that Jy M{v')e^ dv' is uniformly bounded on [0, T] x R" x V, 
we deduce that cp'^ does not depend on v. 

Let -0" e (M+ X R") be a test function such that (p^ — tp^ has a local maximum at (t", a;°). We want 
to show that tp^ is a subsolution of ([6|), yielding that (p^ is a viscosity subsolution The supersolution 
case can be performed similarly. Thereby, we define a corrective term rj not depending on e: = ilf' -\- erj. 
The corrector -q is defined up to an additive constant. We choose the renormalization M{v')e~^ dv' — 1. 
We define 77 as follows, 

V(w,v')eyxy, e''^*'"'") = (v'-w) •V,V^°(t,x). (16) 

The corrector rj is well defined. In fact, we can choose any vq ^ V and define e'''^*'^'"-' = fJ-o + {vo ~ v) ■ 
Vx'4'°{t,x). There is a unique positive fio = e^'-*'^''""^ under the condition Jy M{v')e~^ dv' = 1. 

The uniform convergence ensures that ip^ — has a maximum at {t^ ,x^ ,v^), where (i^,a;^) is close to 
(i°,a;°). As F is a bounded set, the sequence {v^) has an accumulation point, say v* . We can extract a 
subsequence (without relabelling) such that (i^, x*^, w*^) — (i°,a;°,w*). We have at (i^,a;^,u^): 

1 - dtip'' - v'' ■ V^V^ = 1 - dtip^ - v^ ■ V^^^ = / M{v')e''^^dv' . 

Jv 

From the maximum property of (i*^, x^ ,v^), the last inequality implies at this point : 

1 - 5tV'(i', x',v') - ■ V,^/;"(^^ x^ v') > f M(w')e''^*''"''"'^~''^*°'"''"'^ci^;'. 

Jv 

Passing to the limit, we obtain at {t'^,x'^): 

1 _ dtip°{f,x°) - V* ■ Vx^°{f, x°) > [ M(v')e''(*'''^'''"*)-''(*"'^°''''^du' = e^^*"'^"-"*) . 

Jv 

From the very definition of the corrector rj (jl6l) . this writes also: 

^veV, 1 - 5*0° (t", x") - V ■ Vxi^\t\ x°) > e''^*"'^"'") . 
Therefore we obtain at point {t^ ,x^)^ 

dv< I M(«)e-''(*°-°-)rf. = l. 



jy l--dt'ip°{t°,xO)-vWxip°{t°,xO) 
We conclude that tp° is a subsolution of ([6]). 



□ 
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